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SYNTHESIS AND POLYPHASE IMPLEMENTATION OF 
WIDE-BAND LOW-PASS CIRCULAR 2D FIR FILTERS 


Doru Florin CHIPER!, Radu MATE? 


Rezumat. Se propune o tehnica de proiectare analitica in domeniul frecven{a pentru filtre 
2D de tip FIR trece jos de banda larga cu simetrie circulara. Proiectarea este bazata pe un 
fitru prototip trece-jos, ideal de maxim plat, cu banda de trecere specificata, derivat din 
functia tangenta hiperbolica. Prototipul ideal este aproximat cu un polinom trigonometric, 
printr-o schimbare de variabila si dezvoltare in serie Chebyshev, astfel rezultand raspunsul 
in frecventa factorizat al prototipului 1D. Aplicand prototipului o transformare de 
frecventa specifica, derivata din transformarea generala McClellan, va rezulta direct 
raspunsul in frecventa al filtrului 2D circular dorit, de asemenea factorizat. Filtrul 
proiectat are o forma precisa, cu distorsiuni neglijabile. S-a propus de asemenea si o 
implementare eficienta la nivel de sistem pentru filtrul proiectat, utilizand metoda de 
descompunere polifazica si filtrare pe blocuri, ce duce la o structura de filtrare eficienta, 
cu complexitate aritmetica redusa si un grad ridicat de paralelism. 


Abstract. An analytical synthesis technique in the frequency domain is proposed for a 
particular class of 2D filters, namely circular wide-band low-pass FIR filters. The design 
starts from a low-pass prototype filter which is ideally maximally flat, with specified 
bandwidth, based on hyperbolic tangent function. This is approximated as a trigonometric 
polynomial using a change of variable and the Chebyshev series, thus obtaining the 
factored frequency response of the FIR filter prototype. Applying a specific frequency 
mapping derived from the more general McClellan transform, the frequency response of the 
desired circular 2D FIR filter results directly, also in factored form. The designed filter has 
an accurate shape, with negligible distortions. We also proposed a computationally efficient 
implementation at system level, based on polyphase decomposition and block filtering 
approach, which leads to a filtering structure with low aritmetic complexity and a high 
degree of parallelism. 
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Introduction 


Synthesis techniques and implementation structures for two-dimensional filters 
are fundamental research topics in the vast domain of digital signal processing. 
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Some currently used design methods rely on global numerical optimization, while 
others are analytical, employing some frequency transformations applied to 
various 1D filter prototypes, deriving directly the 2D filter with a desired 
frequency response [1]. A specific tool for 2D filter design is the McClellan 
transform [2-4]. A specific class of 2D filters, widely used in image processing, 
are circular filters, of either FIR and IIR type, developed in early papers like [5-7]. 
Some works such as [8, 9] elaborate novel techniques regarding both design and 
implementation of circular filters. As regards implementation aspects, there have 
been proposed various 2D filter structures, both for FIR and IIR versions. The fast 
block implementation of 2D digital FIR filters was proposed in early papers like 
[10, 11]. A high performance 2D parallel block-filtering system for real-time 
applications was presented in [12]. Previous related papers of the authors on the 
design and implementation of some circular filters are [13-15]. 


In this work we propose an analytical synthesis technique for 2D circular FIR 
filters. The method is based on a prototype filter, to which a specific 1D to 2D 
frequency transformation is subsequently applied. We have chosen a maximally- 
flat, smooth function, easy to approximate by a trigonometric polynomial, as 
shown next. Then, an efficient system-level implementation was developed for the 
designed filter, using a polyphase decomposition of a 2D filtering operation with a 
large kernel size. 


Wide-band low-pass zero-phase prototype for 2D circular filters 


As prototype for the desired 2D circular filter, we choose a continuous, smooth 
function with a high steepness, namely the hyperbolic tangent. Specifically, one 
ideal function for a low-pass maximally-flat filter may have the expression: 


H,(@)= 0.5-(tanh(10-(@+2/2))—tanh(10-(@-7/2))) (1) 


Now we look for a trigonometric series expansion of the ideal function H,(@), 
which has to be an approximation as accurate as possible on the frequency range 
[—z, 7]. The most convenient for our purpose is the Chebyshev series expansion, 
because it yields an efficient approximation of a given function, which is uniform 
along the desired interval. The Chebyshev series of a given function on a specified 
interval can be easily found using a symbolic computation software like MAPLE. 
However, we will need rather a trigonometric expansion of H,(@), namely in 
cos(n@) , rather than a polynomial in the frequency variable @. Therefore, prior 
to Chebyshev series calculation, the following variable change will be applied: 


@ =arccos(x/) <> x = cos(a) (2) 


and thus a general, ideal prototype H,(@) can be further written in intermediate 
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variable x as H,(@) =H, (arccos(x/7)) , then the Chebyshev series approximation 


is found as the polynomial expansion in variable x: 


N 
H, (x) = Hp(x)= oy: XH" = Cy FOXES HOG? tit Cy x” (3) 
n=0 
This results indirectly, using the change of variable (2). For the specific prototype 
(1), the following intermediate function in variable x results: 


H (x) =0.5- (tanh (10 (arccos(x/z)+ n/2)) — tanh (10 (arccos (x/z) -x/2))) (4) 
As an odd-parity function, it has the following Chebyshev series approximation, 
in odd powers of x: 
H,,(x) = 0.499726 + 1.086074 - x —0.910301- x? +.0.441843- x° 
—0.113842- x’ +0.015767-x° —0.0011094- x'! +.0.00003111- x! 


The number of terms N=13 is chosen to ensure a desired precision (specified by 
the maximum error, imposed of value 0.06 in our case). 


(5) 


Next, substituting back x= cos(@) , the following frequency response of the low- 
pass prototype filter, in factored form, is derived, where y=cosa: 
H ,(@) =89.741-(y +0.99053)(y + 0.943854)(y + 0.849386)(y + 0.71926) 
(y+ 0.553156) y +.0.369463)(y +0.2044424) .( y? —0.87852- y +0.28381) (6) 
( y? -1.65878- y +0.7454)( y? —2.09279- y+ 1.102997) 


This prototype has the frequency response shown in Figure 1, b and we notice 
it has a relatively small ripple in the passband and stopband. 


H,(@)'?5 ii @y"* 


1 1 


0.8 0.8 


0.6 0.6 


-0.2 -0.2 


3 2 1 0 1 2 3@ 3 2 ‘A 0 1 2 3@ 
Fig. 1. Wide-band low-pass prototype filter 


a. Ideal prototype based on hyperbolic tangent b. Approximated prototype filter 
function 
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Design of circularly-symmetric 2D FIR filters 


As shown before in a particular case, substituting back x= cos(@) in the 
polynomial expression (3), we obtain an approximation Hp(@) of the given ideal 
prototype H,(@), as a trigonometric polynomial in powers of cos@ [15]: 


N 

H,(@) = H,(@)= yh -cos"(@) = by +b, cos @ +b, cos” @+ b, cos’ @+...+ by cos” @ (7) 
n=0 

with by) =cy and b, =¢, ya According to the fundamental theorem of algebra, 


the polynomial (7) can be factored into first and second order polynomials in 
cosa , as follows (where n+ 2m=N , the filter order): 


H(a)=k-| | coso+a,)-| | (cos* o+a,,cos@+a,,) (8) 
i=] 


j=l 


We describe in this section an efficient design technique for 2D circularly- 
symmetric filters, based on the previous 1D filters, considered as prototypes. 
Given a prototype with transfer function H,(@), the 2D circular filter function 


H_(@,,@,) results applying the mapping @ > Jo, + o, : 


H e04,0))= Hp ( Joe +08 | (9) 


The currently-used approximation of the 2D circular cosine function cosa + @; 
is given by the 3x3 array [4]: 
0.125 0.25 0.125 
C=| 0.25 -0.5 0.25 (10) 
0.125 0.25 0.125 


such that we have the following approximation, a particular case of the McClellan 
transform [2-4]: 
cos /@, + @, = C(@,,@,) = —0.5+0.5(cos @, +.cos@,)+0.5cos@,-cos@, (11) 


In order to obtain a filter with circular symmetry from the factored 1D prototype 
function, we simply replace cos@ with the circular cosine function (11) in the 
prototype frequency response (8). Corresponding to this factored expression, the 
filter convolution kernel A results as: 


A=k-M, *M, *...*M, *N, *N, *...*N,, (12) 
where M; (i=1...2) are 3x3 arrays and N; (j=1...m) are 5x5 arrays, given 


by the expressions: 
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M, =C+a,-Ao, (13) 
N, =C#*C+a,,-Cy +a,;- Ay (14) 


where Ao, is a 3x3 zero array and Ay is a 5X5 zero array with the central 
element equal to one; Cy is a 5x5 array obtained by bordering C (3x3) with 
zeros. The above expressions correspond to the factors in expression (8). 

In our case, the component matrices (7 matrices of size 3x3 and 3 matrices of size 
5x5) result according to expressions (13) and (14). For instance, M, is given by: 


0.125 0.25 0.125 
M, =| 0.25 0.49053 0.25 (15) 
0.125 0.25 0.125 


The other 3x3 matrices, M,-M, derive from matrix C given by (10), substituting 


the center element with the values, respectively: 0.443855, 0.349386, 0.21926, 
0.053156, -0.130537, -0.295557. The 5x5 matrices N,, N, and N; are the 


following: 


0.015625 0.0625 0.09375 0.0625 0.015625] 
0.0625 -0.10981 —0.34463 -0.10981 0.0625 

N, =| 0.09375 —0.34463 1.285567 -0.34463 0.09375 (16) 
0.0625 -0.10981 —0.34463 -0.10981 0.0625 
0.015625 0.0625 =0.09375 =: 0.0625 0.015625 | 


0.015625 0.0625 0.09375 0.0625 0.015625] 
0.0625 —-0.20735 -0.53969 —-0.20735 0.0625 

N, =| 0.09375 -0.53969 2.13729 —-0.53969 0.09375 (17) 

0.0625 —-0.20735 -0.53969 —-0.20735 0.0625 

0.015625 0.0625 0.09375 0.0625 0.015625 | 


0.015625 0.0625 0.09375 0.0625 0.015625 
0.0625 -0.26159 —-0.64819 —-0.26159 0.0625 
N; =| 0.09375 -0.64819 2.71189 —0.64819 0.09375 (18) 
0.0625 -0.26159 —0.64819 —0.26159 0.0625 
0.015625 0.0625 0.09375 0.0625 0.015625 


The frequency response and its corresponding contour plot for the designed 
circular wide-band low-pass 2D FIR filter, derived by frequency mapping from 
the prototype in Figure 1(b), with the specified bandwidth, are displayed in Figure 
2. It can be noticed that the filter has a relatively precise shape in the frequency 
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plane and a given ripple in both the pass band and the stop band, corresponding to 
the imposed specifications. 


H.(Q,,@) I 
14 f ] 
O, | | 
0.5 J | 
| 
oJ. & | 
0 ie = — aa a: —iT | 
QO, Pe 0 M, —7T QO, 1 
Fig. 2. Wide-band circular FIR filter 
a. Frequency response of the 2D circular FIR b. Corresponding contour plot in the frequency 
filter plane 


Polyphase implementation of the 2D circular filter 


In this section an efficient implementation of the designed FIR circular filter is 
proposed, based on a polyphase decomposition of a 2D filtering operation with a 
large convolution kernel (27x27). In order to tackle this high complexity problem 
we employ a block processing approach [16] and a polyphase decomposition. 
Using related ideas as in [13], we have first derived a 2D 3x3 filtering algorithm 
that is presented as follows. 

To achieve this task, the kernel of the 2D filter obtained from design and the input 
image to be filtered are decimated by factors 3 and 9, respectively; after this, a 
polyphase filtering technique is used. Through this approach, three partial output 
component images are calculated in parallel, namely Y,, Y, and Y, given by 


equations (19), (21) and (23): 


Ao Ao T 
A’ O,, O A, O,; -A; O,; O 
Y)=|035 Ag | ae alia et cae “|e (19) 


O35 Ao Ay Ay Ag O65 Ar Ar O65 O65 
where the vector Vo is: 
Vy=[1 1/2 1/2 11 1 1/2 1/4 1/4 1/2 1/2 1/2]7 (20) 
Similarly, we have: 


Ay 03, 


Al =Al At O,; A, -A,; O.; O 
Y,=|-Ay Ag |-diag(V,)xdiag : - . nt ‘ ae ; : oe “|i (21) 
A, O 063 Ap O63 Os Ar Ary O65 O65 
3,6 
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Yea [2? VA Tay To ASA A. Yh a? (22) 
Y, = re - shade | Ao alee Ay Oss Ai 68 923) 
Sac Ay O63 O63 Ao Ocs -A; -A, A, Ay 
V, =[1 1/2 1/2 11141 1/2 1/2 11 |’ (24) 
In the above matrix relations, the 3x6 and 6x5 block matrices are: 
1 0-1 00 
0 1 1 00 
Lied Det FOr O Gy 2g? EOF 6 
A,=|0 1 -1 1 1 0 hae | ee ae ee er (25) 
01 10-11 a eee ee 
0 -1l -l 1 1 


and O;,, O.;,0,,5 are zero matrices of size 3x6, 6x3, and 6x5, respectively. 
Summing up the obtained partial results Y), Y, and Y,, the following output 
vector Y containing 25 samples of the filtered image is obtained: 
T 
Y=Y+Y,+Y, = [Yoo Yr Yoo No YN No Yo Yo Yoq | (26) 
Also, we have: 


h=[Noo My, Mon Myo yy M2 My My Ayg | (27) 


X2p = [00 Xo1 02 %03) X04 Ato M33 X34 X4g X41 X42 X43 x44] (28) 
The above algorithm for a 2D FIR filtering operation has been derived using a 
subexpression sharing technique in order to reduce the number of operations due 
to the fact that there are redundant operations in the direct implementation using 
2D convolution. Thus, in direct 2D convolutions there are input data blocks that 
are overlapping, and avoiding these redundant operations we can obtain important 
savings in the number of arithmetic operations. Using a subexpression sharing 
technique, the above algorithm resulted for the 2D FIR filtering operation. 
In order to extend 2D filtering operations to a 27x27 kernel we will use next a 
polyphase decomposition with a decimation factor of 3 and thus we obtain 9 
submatrices. For the input matrix we have used a decimation factor of 5 and we 
have obtained 25 submatrices for a 45x45 input matrix. 
The equation (27) can now be reformulated for a polyphase implementation by 
replacing the vector h with the vector H defined as below: 


T 
H =| Ho Ay, Ay Aig A, Ay Ay) Ay, Hy | (29) 


The vectors Hoy,..., 42 for the size 27x27 of the kernel matrix are defined as: 


0 1 2 3 4 oD) 6 7 8 
Hy =| Hoo Ay A Ay A Ay A Ay Ayo | (30) 
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where the components of vector Hog are given in the group of equations (31): 
Hoo =[o0 'o3 Moe oo Mor tos Mors Moo Mor | 
Hy =['y0 Mya Mae Meo Marr Mais Mars Mor Maza | 
He ae [Neo he 3 Neo Neo Ne12 Nes Neg No.1 hes | 
Hoo 7 [he o Myx gg Mg Mgr Mois Moig Mo,21 hy 04 | 
Hoo = ls hyn 3 hi.6 Iyo.9 Iyo32 Myo15 hyo18 Iyo>1 Ino 04 | (31) 


5S 
Ay = 
6 _ 
Ay = 


Ai 


8 _ 
Ay = 


his, hs his,6 


higo hig3 hig. 


Myo M3 Mig 


hy4.0 Iy43 hy4 6 


his. Iys,12 Isis hiss Iy5,>1 his | 


hig Iygio Migs ligig Iyg.o1 hiss | 


hy19 hy) 12 Iy\ 15 yi 18 hy ho 2 | 
In4.9 Iy4 19 Iy4 15 Iy4 18 Iy4 91 okay | 


Similarly, the other vectors and their components are displayed in the following 


groups of expressions: 


0 1 2 3 4 5 6 7 8 
Ao =| Aoi Ay, Ao, Ao, Aor Hor Ao Aoi Ay, | 
where 
oie 
Ho =| Mo. M4 07 Mo10 "13 Mo16 0,19 %,22 hy. | 
a 
Ao, =|'31 'a4 Ia Mao Ia M16 Ib19 Ia, hs. | 
io 
Ao, =| N61 Nea he7 N6,10 he.33 N6.16 N6,19 Ne22 hes | 
aaa 
Ho =| M1 M4 Mo7 Mo19 Mo,13, Mo,16  Mo19 No.2 hy os | 
ee 
Ay, =| M21 Nya Inna Iyn10 Myo 43 Ii16 Ino19 Nyy 29 hyo s | 
ee 
Ao, =| 51 isa Iisa is10 ns.13 Nis.16 i519 Ns,09 hss | 
oe 
Ao, =| ns ga liga igo igi3 lige Nigig Nigoo figs 
eae 
Ay, =| oi Nyi4 Iyy7 y110 Nyi33 hyi16 Ny119 Ny1.20 hors | 
ee 
Ao, =| Moa Nog 4 Nyy Ny4.10 Iya3 Ny4.16 Ny4.19 Ny4.99 igen 
also 
(0) 1 2 3 4 5 6 we 8 
Hy | Ao ahS) A) Ap A Ay Ao) Ao Ay, | 
where: 
0 
Hy | las liga: Tiles May: Naglae’ Dense Magee Tas. “Toe:| 


Ins Fag Ion Ia lar an Mas a2. | 


(32) 


(33) 


(34) 
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aT 
Ao =| Neo Nes 6g Neotr Nera Moz "6,20 "6,23 her | 
36 
Hop =| No.2 ys Nog Mgiy Moja Moi7 M290 Io,23 hy a6 | 
4 
Ae =[A22 Mins Mog Mo Mora M217 M220 "h2,23 hyo 06 | (35) 
. oT 
Hop =| 5,2 iss Msg Msi Msia Msi7 M520 5,23 his | 
6 oT 
Ao, =| Ihg.r igs Mgg Meir Msia Msi7 M20 823 hisze| 
7 oT 
Ao =| Mai2 Ays Mg Mia Moa Misr Min Moi hoi | 
sor 
Aon =| Ina Mygs Mogg Moan, Moaia Moan Iha20  ha3 hogze | 
also 
[ 770 1 2 3 4 5 6 7 8 
Ay =| Ayo Ayy Hy Hig Hy Hig Aion Ato Ay | 
where 
oT 
Ay =| ho ns Me Mo Mar Mas Mas M21 I>, | 
cor 
Ay =| hao Ig3 Mg Mag Mair Mais Marg Mari hao | 
+ oT 
Ay =| ho ny Mne Myo Mar Mais Mais M21 Ino | 
3 oT 
Hin =| noo Ios Moe Mos Moir Mors Mors “4o21 Ino | 
4 
Ai =| hia 33 M36 M39 Maire Mais Misig 3,21 Iya | (36) 
sor 
Aig =| he. ies yee Moo Moir Mois Mois "he21 Iie | 
6 oT 
Ay =| hoo Nos Moe Moo Moir Moras Mors /ho21 hiy.o4 | 
7 oT 
Ay =| hy0 Myx M6 M9 Mir Mis M18 M221 hyo 24 | 
er 
Ay =| Nos. Iy53 My56 M59 M532 ysis Mosag (5,21 Ios. 04 | 
and 
[ 770 1 2 3 4 5 6 7 8 
A, =| Ai A, Ay Ay Ay Ay, Ay, Ay Ay, | 
where 
— 
A, =| hy ha M7 Mio Mas Mie Mao M2» Is | 
Lor 
A, =| hay Mgg Igy Mgin Mai3 Maio Mano Ma20 Iyos | 
> oT 
A =|hy1 na yy Inyo ys Mie Mio 4x Ins | 
+ oT 
Ai, =| Mo. tions Moz Moro “013 M016 Mo19 0,22 Ios | 
4 
Hy =| hay liza M37 Maio Mais aie M319 3,22 Iyas | (37) 


=| Mn6.1 Ihe. hie hy6,10 Ih.13 Ii6,16 Ii6,19 Iyo.22 Ios | 
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6, 
Hie tis Ino 4 Ino7 Iyo10 Iyo13 Iho16 Iyo19 Iyo,29 hos | 


a 
Hi, =| 2,1 hyp 4 hy 7 hy 19 hy 13 hy 16 hy 19 hyp 99 hyp.9s | 


Hy, = Mos, Mos Mos7 M510 In513, Ia5,16 Mas.i9'n5,22 hys.95 | 

and 

Fy =| Ap Hip Hyp Hi, Hi H}, Hi Hi Hy | (38) 
where 


oO _ 

Ay = hy» hs hg hyiy Iya hyi7 hy 0 hy 93 Ie | 
eee 

Hy) =| hy» hy hy g Igy hgi4 hg7 hg 00 hg93 hye | 
27 

Ay) =|lh2 hy 5 hy Iya) My ig hy 17 hy 0 hy 93 In 36 | 


Bn 
Ay, = Ino2 hos hog yo Ino,14 Iyo.17 hio,20 Iio,23 hio,26 


4 
Ais =| hg. Iy3,5 hy. y3iy My314 My317 M1320 My3.03 Iya. | (39) 


sf 
Ay, =| hhe.2 hie.s hieg Mie. Iy6,14 Iy6,17 Ii6,20 Ih6,23 I16,26 


6. 
Ay = Io Ino Ihog Myo Iyo14 Iyo17 Iyo,20 Io 23 Iho 26 


7 
Ay, = hyp > hy 5 hy g hyp 1] hy 14 hyp 17 hyp 09 hyp 93 hyo 26 | 


8 
Ay = | 95,2 Nhs 5 Iys 8 ys Nys14 Nys17 Nys5,09 Iys5,93 hoes | 
and also 
(0) 1 2 3 4 5 6 7 8 
Ay =| Hoy Ay A459 5 5 59 5 A459 Hy | (40) 


where: 


Hin =| tg In3 Ing Ing Injr Inis Mors Mor hy o4 | 
Hy =| Aso Mis3 M56 Its9 M512 M515 Ms1gs,21 hs. | 
Hi | fee. Pia te Tay. Teas Tate ag ge ha, 
Hx =| M10 Mis Mio Mio Mure Mis Maas M121 Ina | 


4 
59 a [4140 ha3 Iha6 Iya9 yao Maas Iyaig Muar Ingo | (41) 


5 _[ 
A) =| In79 73 hy7,6 hi7z9 hyz12 Iyz15 hy718 Iy7,21 Ing | 


6 
Ay = hyo hyo3 hyo6 hyo.9 Myo hyo15 hyo18 hyp 91 fa5i| 


4 = 
Ay =| M3. hy3.3 hy3.6 hy39 hy319 hy315 hy3 18 hy3.91 Iy3.24 | 


8 
Aq =| 46.0 Iye3 hy6.6 Iy69 Iy6.19 hy6.15 Iy6.18 hy6,>1 Iyg.24 | 


0 1 2 3 4 5 6 7 8 
Hy, =| Ha, Hy, Hy, Hx, Hy, Hy, Hy, Ay, H} | (42) 
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0 _ 

Ay, =|Iyn, Ing 
Les. 

Ay, =|h5,  Is,4 


2 
Ay, = he hg 4 


Ing Pos Ioas Ios Maas Inze Ia | 
hs 4 hs 19 hs13 hs 16 hs 19 hs 99 hss | 


ay Pero Meas Mets Meas Maar Meas | 


Backs 
A, =| In Ini 4 fyi Ay 1.10 fyi33 My 116 My 119 hy 1,29 hiias | 


dite 
Ay, — May Iga 
H>,=|ho, h 

21 17,1 17,4 
HS, = 

a1 = Nyo1 Nyo.4 
Hi =|h 

21 23,1 Iy3.4 


So. 
Ay, = | 261 hye, 4 


0 1 
Ay, =| Hy Ay 


oT 
Ay, = hy» hy 5 


I 
Ay =| hs hs.5 


2 
Ay, = hg» hg 5 


Ia7 M40 Iya3 Iha.16 Ia19 I4.o9 Ings | (43) 
hy,7 hy7,10 hy713 hy7,16 hy719 hy7,29 hy7,95 ] 

hy7 hyo10 hyo13 hyo16 hyo19 hyo,99 hyp 25 | 

hy37 hy310 hy313 hy3.16 hy319 hy3,99 hy3 95 ] 

hyo7 hy610 hy6.13 hy616 hy6.19 hy6,99 hogs | 


2 3 4 5 6 7 8 
Any Ay Ay, Ay, Ay, Any Hi, | (44) 


Ing Poa toa Iaay Po Ios Ih p6 | 
hs g hs 44 hs 14 hs17 hs 90 hs 93 hs. | 
fg Morr Pera Marr Pern Pens Pore | 


3 
Ay =| hy hy 15 hy Ayia Iya Myia7 hy 1,20 hy 1,23 hiss! 


4 _ 
Ay) — [ thas has 
Heal he: Hh 

1) a 17,2 17,5 
HS, = h 

22 — Nyo.2 20,5 
Hi,=|h 

22 — L 23,2 Iy35 


Boe 
Hy = hyo hye, 5 


hag Mya Maia Muay Iy4,20 Iya.o3 Ina | (45) 
hz Iya Iz 14 Iyza7 hy7,20 hy7,33 hia 26 | 

hyo hyo 14 hyo14 hyo17 hyo29 hyo3 Iy9.26 | 

hy3.g hy31) hy314 hy317 hy3 99 hy3.3 In3.36 | 

lye hy614 Iy614 hy617 hy6,30 hy6.93 Iy6.26 | 


For the input matrix of size 45x45 we have organized the data as shown below. 
We also replace the vector x,, with the vector X,, defined as follows: 


T 
Xap =[Xoo Xo Xo Xo3 Xo4 X19 Xi oS X33 X34 X49 Xa X49 X43 Xu] (46) 
The vectors Xo), X01... 44 are defined below: 
8 7 6 5 4 3 2 1 0 
Xo = [ Xt Xoo Xoo X00 Xoo X00 X00 Xoo Xo | (47) 


where X(),.-,X) are derived from the following vectors Xai ee Xavo . whose 
elements are taken in reversed order: 
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oT 
Xd = 


Xapo = 
Xabo = 
ae = 
Xavo = 
xa = 
Xaby = 
Xajy = 


8 _ 
Xd = 


X00 %o5 *o10 *015 0,20 0,25 0,30 0,35 ib | 


X50 6%5,5 510 %515 5,20 5,25 %530 %5,35 5,40 | 


[*10,.0 ios *10,10 1015 10,20 10,25 10,30 10,35 4646 | 


15.0 15,5 15,10 1515 15,20 15,25 15,30 715,35 %5,40 | 


| *20.0 20,5 20,10 2015 20,20 20,25 *20,30 ~—*20,35 en 


%75.0 25,5 2510 2515 25,20 25,25 25,30 25,35 5,40 | 


| *30.0 %30,5 30,10 30,15 30,20 *30,25  *30,30 30,35 545 | 


%35,0 %35,5 35,10 3515 35,20 35,25 35,30 35,35 -*35,40 


| *40.0 40,5 40,10 4015 40,20 40,25 40,30 %40,35 en 


Similarly, we have: 


Xo. = 
with 


an. 
Xd, = 


Xap = 
Xap) = 
Xap = 
Xa, = 
Xap) = 
Xap) = 
Xagy = 
Xap) = 


and so 


8 7 6 5 4 3 2 1 0 
[ X61 Xo Xo Xo Xo1 Xo1 Xo Xo1 Xi | 


%o1 %0,6 *011 016 0,21 0,26 0,31 0,36 aia | 


| %5.1 %56 %511 516 75,21 5,26 5,31 5,36 5,41 | 

Ear Moo 1011 *io16 10,21 *10,26 ©=*10,31 10,36 Moai 
En M56 1511 71516 15,21 15,26 %15,31 15,36 “15,41 | 
Ea %20,6 20,11 2016 20,21 *20,26 20,31 20,36 *20,41 ] 
En %25,6 *2511 *2516 25,21 25,26 25,31 25,36 © *25,41 ] 
En %30,6 *30,11 30,16 30,21 30,26 30,31 *30,36 © 30,41 | 


351 %35,6 3511 35,16 %35,21 *%35,26 ©%35,31 35,36 ©=-*35,41 


| *40,1 40,6 40,11 40,16 40,21 40,26 40,31 *40,36 | 
on, until the last vector: 

[ y8 7 6 5 4 3 2 1 0 

| X44 X4g Xgy Xyyg Xp Xe Xe Xyy x4, | 


X44 X49 %414 419 %424 %4,29 %4,34 4,30 aie | 


|%9.4 *99 *o14 %o19 9,24 %9,29 7934 9,39 X44 | 


=|[%144 7149 A414 A419 714.24 A1429 M1434 M1439 4.4 | 


=|%19,4 %19,9 1914 *1919 719,24 719,29 %19,34 19,39 io.a4 | 


(48) 


(49) 


(50) 


(51) 
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‘cr 
Xayy =|%24.4 %249 %2414 %24,19 *%24,24 24,20 24,34 24,30 24,44 


(52) 


Gi cs 
Xay4 =| %344 %34.9 3414 %3419 %34,24 %X34,29 34,34 %X34,30 34,44 


5 7 
Xay4 =| X29,4 %999 2914 2919 29,24 29,29 29,34 29,39 X20,44 | 


vine 
X44 =| %X39,4 %309 30,14 3019 %30,24 © %30,29 30,34 30,30 © 30,44 


Xai =|%444 %449 Xaa14 %X4a19 %44.24 X44.290 X44.34  X44,30 is de 

Due to the lack of space, we have shown only the first two elements and the last 
one. Applying the decimation with 5 instead of the input matrix of 45x45 we 
obtain 25 matrices of size 9x9 that have been reformulated as 25 vectors with 81 
elements. Thus, we have decomposed a 2D FIR filtering operation with a 27x27 
kernel and a 45x45 input matrix into 25 inner product operations. 


Conclusions 


The analytical design procedure described above is simple and efficient, leading 
to 2D filters of accurate shape. The prototype chosen here is a wide-band low-pass 
filter with a specified bandwidth, to which a simple frequency mapping is apllied, 
yielding directly the factored frequency response of the desired 2D circular filter. 
The method is more general and can be also applied to obtain 2D filters with other 
shapes. Also, an efficient implementation structure at system level for 2D FIR 
filters has been proposed. Employing the polyphase decomposition and the block 
filtering approach, a given 2D filtering operation of an image with a large size 
filter kernel can be achieved equivalently as a sum of filtering operations with 
several smaller size kernels. Instead of the redundant and resource consuming 
direct convolution of an image with a large kernel (in our case, 27x27), the 
decimation by given factors and block filtering are used, such that computational 
complexity of the overall filtering task decreases significantly. 

In further work, we will focus on determining the optimal decimation factors for 
the filter kernel and input image in order to obtain the most efficient structure. 
Also, we envisage developing such implementation structures for other types of 
2D filters, all having a low arithmetic complexity and high parallelism. 
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